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ABSTRACT 



Aims. To study the variation of the angular momentum and the rotational kinetic energy of the Sun, and associated variations in the 
gravitational multipole moments, on a timescale of the solar cycle. 

Methods. Inverting helioseismic rotational splitting data obtained by the Global Oscillation Network Group and by the Michelson 
Doppler Imager on the Solar and Heliospheric Observatory. 

Results. The temporal variation in angular momentum and kinetic energy at high latitudes (> tt/4) through the convection zone is 
positively correlated with solar activity, whereas at low latitudes it is anticorrelated, except for the top 10% by radius where both are 
correlated positively. 

Conclusions. The helioseismic data imply significant temporal variation in the angular momentum and the rotational kinetic energy, 
and in the gravitational multipole moments. The properties of that variation will help constrain dynamical theories of the solar cycle. 

Key words. Sun: interior - Sun: rotation - Sun: activity - Sun: helioseismology 
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1. Introduction 

Helioseismic data from the Global Oscillation Network Group 
(GONG) and the Michelson Doppler Imager (MDI) have been 
used to infer the rotation rate in the solar interior (Thompson 
et al. 1996; Schou et al. 1998). With the accumulation of these 
data for about 11 years it has now become possible to study 
its temporal variations. Previous studies have revealed a dis- 
tinct pattern of bands of faster and slower than average rotation 
(Howe et al. 2000a; Antia & Basu 2000), which at low latitudes 
move towards the equator, and at high latitudes move towards 
the poles (Antia & Basu 2001). This pattern, which is similar to 
the torsional oscillations observed at the solar surface (Howard 
& LaBonte 1980), has been well studied, and it has been found 
to penetrate through much of the body of the convection zone 
(Vorontsov et al. 2002; Basu & Antia 2003; Howe et al. 2005, 
2006). Below the convection zone, however, there is no generally 
accepted unambiguous detection of significant temporal varia- 
tion in the rotation, although Howe et al. (2000b, 2007) have re- 
ported a periodic variation near the equator between r = 0.65/? Q 
and r = 0.75/? Q in the early years of the last sunspot cycle with a 
period of 1.3 years, and Gough (2007) has reported that there is 
evidence that the oscillation penetrates even more deeply. Other 
studies (Antia & Basu 2000; Basu & Antia 2003) have not con- 
firmed such a variation. 

A direct consequence of temporal variation in the angular 
velocity is that global quantities like angular momentum and ki- 
netic energy of rotation must also exhibit variation with time 
(cf., Komm et al. 2003). Because of the steep increase in den- 
sity with depth, the major contribution to these global quanti- 
ties is naturally likely to arise from the deeper layers, where a 
temporal variation in rotation has not been determined reliably; 



Send offprint requests to: H. M. Antia, e-mail: antia@tifr . res . in 



and as a result, it is difficult to draw any firm conclusion about 
the global variation of angular momentum and kinetic energy. 
However, following Komm et al. (2003), we may consider the 
contributions from different layers within the convection zone 
separately. We hope that such a study will shed some light on 
the rotational dynamics and its role in driving the activity cy- 
cle. Indeed, the study we present here is very similar to that of 
Komm et al. (2003), although we are now able to extend it over 
an entire solar cycle. Apart from angular momentum and kinetic 
energy, we study also the temporal variation of the gravitational 
quadrupole and higher-order multipole moments due to the dis- 
tortion in the solar figure induced by the known rotation. 

In this work we concentrate principally on temporal varia- 
tions in angular momentum and rotational kinetic energy in dif- 
ferent regions in the solar convection zone, and compare their 
variation with the solar activity. We describe in Section 2 the 
basic technique and the data used in the study; the main results 
are described in Section 3. In Section 4 we present the variations 
in the multipole moments of the external gravitational potential. 
We draw our conclusions in Section 5. 



2. The data and their analysis 

We use the data from GONG (Hill et al. 1996) and SOI/MDI 
(Schou 1999). Each data set consists of mean frequencies of 
different («, /) multiplets, and the corresponding splitting co- 
efficients. We use 110 temporally overlapping data sets from 
GONG, each covering a period of 108 days, starting from 1995 
May 7 and ending on 2006 May 20, each set being displaced 
by 36 days from its predecessor. The MDI data consist of 49 
contiguous data sets, each covering a period of 72 days, start- 
ing from 1996 May 1 and ending on 2006 May 20. These data 
cover the rising phase of the cycle 23 and most of the descending 
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phase as well. We use a 2D Regularized Least-Squares (RLS) in- 
version technique in the manner adopted by Antia et al. (1998) 
to infer the angular velocity in the solar interior from each of the 
available data sets. 

To study the temporal variation in the angular velocity we 
look at the residuals obtained by subtracting from the angular 
velocity at any given time its temporal mean Qo(r, 0) 



6Q.(r, 0, t) = Q(r, 0, t) - Q (r, 9) 



(1) 



with respect to the spherical polar coordinates (r, 0), where t is 
time. The averaging to obtain Qo is over the entire time inter- 
val in which the data are available. Both the temporal mean and 
the residuals were calculated separately for the GONG and MDI 
data. The residuals, which essentially define the temporally vary- 
ing component of the angular velocity, were used (together with 
the means) to calculate the variations with epoch in the kinetic 
energy and angular momentum in different regions of the Sun. 
For example, the change in angular momentum of the Sun be- 
tween spherical surfaces at radii r\ and r 2 can be written as 

6J(t ; ri , ri) = 2n f pr 4 dr f SCl(r, 9, t) sin 3 d0 , (2) 

Jr, Jo 

where p(r) is the density in the Sun. The integrand in the outer 
integral defines the contribution to the angular momentum per 
unit radius from an infinitesimal spherical shell around radius r, 
namely dSJ/dr :- d5J/dr 2 \ ri=r . Similarly, we can also calculate 
the temporal variation in the kinetic energy of rotation from the 
expression 

ST{t;n,r 2 )^ 

pr 4 dr I Q (r, 0) (5Q(r, 0, t) sin 3 d0. (3) 

r, JO 

In evaluating expressions (2) and (3) we have used the density 
p(r) from a standard solar model. 

In addition to integrating over the entire latitude range we 
can consider also contributions from different latitude intervals. 
The observed variation in angular velocity exhibits different be- 
haviour in the low- and high-latitude regions. At low latitudes 
the bands of more rapidly and more slowly rotating fluid mi- 
grate towards the equator, while at high latitudes the bands move 
polewards (Antia & Basu 2001; Vorontsov et al. 2002; Howe et 
al. 2005; Basu & Antia 2006). We therefore consider the contri- 
bution to the integral separately over the low latitudes (< n/4) 
and high latitudes (> n/4) to see how they compare with the 
corresponding global quantities. 

The integrands in equations (2) and (3) are essentially the 
same, save for the factor Qo in equation (3). Therefore because 
the angular variation of Qo is substantial only in the polar regions 
where the factor sin 3 6> is quite small, the temporal variations 6 J 
and 5T are very similar. Consequently in most of the following 
discussion we consider explicitly only the temporal variation in 
the kinetic energy. We compare it with a solar activity index, for 
which we use the radio flux at 10.7 cm as a proxy. 

Apart from kinetic energy and angular momentum, it is 
possible to compute the gravitational quadrupole and higher- 
order multipole moments of the Sun resulting from the centrifu- 
gal force (e.g. Schwarzschild 1947; Sweet 1950; Gough 1981, 
1982; Ulrich & Hawkins 1981; Pijpers 1998; Antia et al. 2000; 
Roxburgh 2001 ; Mecheri et al. 2004). The gravitational potential 
<3>(r, 0) outside the Sun can be written as 
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J 2 k(t)P2k(cOS0) 



(4) 
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Fig. 1. Temporal variation in rotational kinetic energy through- 
out the convection zone. The continuous lines join values of the 
total kinetic energy (in units of the temporal mean) in the convec- 
tion zone at different epochs; the dashed and dotted lines respec- 
tively join corresponding contributions from the low (< n/4) and 
high (> n/4) latitude regions. For clarity, only a few represen- 
tative error bars are shown. The values in the upper panel were 
obtained from the MDI data, the values in the lower panel from 
the GONG data. 



where J 2 /< are the dimensionless multipole moments and P 2 k are 
Legendre polynomials of degree 2k. We attempt to study the 
temporal variations in these quantities. 



3. Results 

We show in Figure 1 the temporal variation of the rotational 
kinetic energy of the entire convection zone, together with the 
contributions from the high- and low-latitude regions. The up- 
per panel was obtained from the MDI data, the lower panel from 
GONG. An oscillatory variation with about an 11-year period 
is evident, although there is some difference between the values 
inferred from the MDI and the GONG data. 

Figure 2 depicts the temporal variation in rotational kinetic 
energy in different regions of the Sun obtained using the GONG 
data. The magnitude of the absolute variation increases with 
depth, largely on account of increasing mass density. In order 
to get a more uniform comparison amongst different depths, 
the variation relative to the mean value of the kinetic energy 
in each layer is plotted; then the amplitudes are rather similar. 
The contributions from low latitudes (0 > n/4) and high lati- 
tudes (0 < n/4) are also plotted. The relative magnitude of these 
variations is also taken with respect to the mean kinetic energy 
in the entire latitude range. In all cases these variations are sim- 
ilar in magnitude. We note that the variation in angular veloc- 
ity is substantially greater at high latitudes than it is near the 
equator. The effect of that on the variation of the rotational ki- 
netic energy and angular momentum is partly compensated by 
the fact that the high-latitude regions are closer to the rotation 
axis, and contribute less to the moment of inertia, a property 
which is reflected by the factor sin 3 in the integrands in equa- 
tions (2) and (3), which on its own integrates to 5/6 V2 w 0.589 
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Fig. 2. Relative variation with time of the rotational kinetic en- 
ergy in different regions (indicated in each panel) of the Sun, ob- 
tained from the GONG data. The continuous lines indicate the 
variation over entire spherical shells, and the dashed and dotted 
lines respectively indicate the variation in the low- (< n/4) and 
high-latitude (> n/4) regions. For clarity, only a few representa- 
tive error bars are shown. The lowest panel shows the variation 
in solar activity as measured by the 10.7 cm radio flux, over the 
same time interval. 

and 2/3-5/6V2 = 0.077 in the low- and high-latitude regions 
respectively. 

At the bottom of Figure 2 is plotted the 10.7 cm radio flux, 
which is an indicator of solar activity. It is evident that in the 
outer layers (r ^ 0.9R Q ) of the Sun the kinetic energy is greatest 
at epochs of greatest activity, both in the low-latitude and in the 
high-latitude regions. At intermediate depths (0.7 55 r/R S5 0.9) 
the high-latitude correlation with activity persists, but at low lat- 
itudes the correlation is reversed. This trend is seen more clearly 
in Figure 3, in which the correlation coefficients are plotted as 
a function of radius, both for entire spherical shells and for the 
high- and low-latitude regions. The amplitude of the variation is 
the greater in the low-latitude region, so the behaviour there is 
an indicator of the behaviour of the kinetic energy, and the an- 
gular momentum, over the entire spherical shell. At the greatest 
depths it is more difficult to see the trends in Figure 2 because 
the uncertainties in the inversions are comparable with the trends 
themselves. However, the covariance of the kinetic energy with 
the activity formally changes sign at r O.7R , at both low 




Fig. 3. Correlation coefficients between the temporal variation 
in kinetic energy of rotation obtained using GONG data and the 
10.7 cm radio flux, plotted as a function of radius. The contin- 
uous curve is the coefficient for an entire spherical shell; the 
dashed and dotted curves are respectively the results for the low- 
latitude (< n/4) and high-latitude (> n/4) regions. 



and high latitudes, beneath which the low-latitude kinetic energy 
variation is again positively correlated with activity and the high- 
latitude variation is negatively correlated. Again, this behaviour 
is visible more clearly in Figure 3. 

These results are consistent with those of Komm et 
al. (2003), who considered the angular-momentum variations 
over entire spherical shells. The relative variations in angular 
momentum and rotational kinetic energy, and also, of course, the 
angular velocity, are of order 0.1%. Given that the total angular 
momentum and rotational kinetic energy of the convection zone 
are J c =s 2.8 x 10 47 gcm 2 s~' and T c =s 4.0 x 10 41 erg, this implies 
that the amplitudes of their variations are about 1 x 10 44 g cm 2 s _1 
and 3 x 10 38 erg respectively. 

The change in sign at r - O.9R of the correlation of an- 
gular momentum with activity at low latitudes could be associ- 
ated with the upward migration of the zonal flow pattern reported 
by Basu & Antia (2003): because the migration timescale from 
intermediate depths is about half a cycle period, the flow pat- 
tern at intermediate depths is anticorrelated with that at the sur- 
face. At high latitudes, the angular-momentum variation is posi- 
tively correlated with surface activity throughout the convection 
zone. Beneath the convection zone the correlation appears to be 
weaker, possibly as a result of the greater uncertainties in the 
inferred angular velocity, and may not be significant. 

To calculate the amplitude of the solar-cycle variation in ki- 
netic energy, we fit a sinusoid with a period of 1 1 years to the 
observed kinetic energy. Since the power spectra at most depths 
show a peak at frequencies corresponding to approximately this 
value, we may be justified in looking for such variation. The 
resulting amplitude of the relative variation is plotted against ra- 
dius in Figure 4. The dips in the curves occur near the depths 
at which the correlation coefficient changes sign (see Figure 3). 
At these depths the amplitude might be expected to be small. 
Beneath the convection zone the rise in amplitude with depth 
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0.0015 




Fig. 4. Amplitude of relative variation in rotational kinetic en- 
ergy, plotted against radius, obtained by fitting a sinusoid with 
a period of 1 1 years to the inferences from GONG data. The 
continuous curve is for the entire spherical shell, the dashed and 
dotted curves respectively are the results for the low (< n/4) and 
high (> n/4) latitude regions. 
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Fig. 5. Correlation coefficient between the the 10.7 cm radio 
flux and the kinetic energy and angular momentum, integrated 
over entire spherical shells, plotted against radius. The solid and 
long-dashed curves are kinetic energy and angular momentum 
inferred from GONG data, the dotted and short-dashed curves 
are from MDI data. 



might be merely a reflection of the increasing errors in the in- 
ferred angular velocity, and is therefore not to be trusted; this 
is suggested partly by the low correlation coefficients evident in 
Figure 3. 

We conclude this discussion by once again comparing infer- 
ences from GONG and MDI. Figure 5 depicts the correlation 
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Fig. 6. Power spectrum of temporal variation in rotational kinetic 
energy in different depth ranges (as marked in each panel). The 
solid lines are for the entire spherical shell, the dashed and dot- 
ted lines are for low- (< n/4) and high- (> n/4) latitude regions, 
respectively. The seismic data were obtained from GONG; the 
results from MDI are similar, and have been omitted for conci- 



coefficients between the 10.7 cm radio flux and the spherically 
averaged rotational kinetic energy and angular momentum ob- 
tained from the two data sets. It is evident that there is very little 
difference between the coefficients for rotational kinetic energy 
and angular momentum. There is some discrepancy between the 
inferences from GONG and MDI in the lower part of the convec- 
tion zone, a discrepancy similar to one that has been found also 
in other works (e.g., Schou et al. 2002). The reason for that dis- 
crepancy is not fully understood. Nevertheless, it does not affect 
the principal conclusions of this paper. 

Howe et al. (2000b) and Komm et al. (2003) have reported a 
1.3-year periodicity in the temporal variation of the angular ve- 
locity at low latitudes in the vicinity of the tachocline. Here we 
examine this finding further, by seeking such a periodicity in the 
angular momentum or kinetic energy. We calculate the discrete 
Fourier transforms of the results illustrated in Figure 2 to see 
whether there is a significant peak in the power spectra around 
a frequency of 0.75 yr _1 . The outcome is presented in Figure 
6. Evidently, no such peak is apparent. The power spectra in the 
vicinity of the tachocline appear to contain no particularly signif- 
icant peak, although those in the upper layers of the convection 
zone do contain an obvious peak in the lowest frequency bin al- 
lowed in the discrete Fourier transform: this corresponds to the 
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Fig. 7. Temporal variation in rotational kinetic energy and an- 
gular momentum in the low-latitude region with 0.70/? o < r < 
0.74/? Q . The continuous lines are the result of using GONG data, 
the dotted lines are from using MDI data. For clarity, only a few 
representative error bars are shown. 

solar-cycle variation. We have looked at many different depth 
and latitude ranges; the results are all similar to those shown 
in the figure, irrespective of whether the seismic data employed 
were from GONG or from MDI. Therefore, apart from a possi- 
ble periodicity characteristic of solar-cycle variation, we find no 
evidence for periodicity at any depth. For illustrative purposes, 
we present in Figure 7 the temporal variation in the low-latitude 
region of the tachocline (Q.1QR Q < r < 0.747? o ) in both ro- 
tational kinetic energy and angular momentum, inferred from 
MDI as well as GONG data. This is the region in which Howe et 
al. (2000b) and Komm et al. (2003) found the greatest amplitude 
of the 1.3-year oscillations. It is clear from the figure that there 
is reasonable agreement between GONG and MDI data, and that 
the variation in angular momentum is similar to that in the rota- 
tional kinetic energy. However, no significant temporal variation 
is evident; in particular, there is no very clear sign of a 1.3-year 
oscillation. 

4. Gravitational Multipole Moments due to Internal 
Rotation 

From our inferences of the angular velocity throughout the solar 
interior, it is possible to calculate the quadrupole and higher- 
order multipole moments of the Sun's gravitational field, and 
study their possible temporal variation. Figure 8 illustrates the 
even-order moments - J12 inferred from both GONG and 
MDI data. To ascertain whether temporal variation in these quan- 
tities is related to solar activity we have calculated the coeffi- 
cients of their correlation with the 10.7 cm radio flux. We iden- 
tify the phase of any possible solar-cycle variation by fitting to 
the moments sinusoids with a period of 1 1 years. The results are 
summarized in Table 1 , in which are listed the temporal averages 
of the moments as well as the relative amplitudes and phases (in 
years, relative to the phase of the 10.7 cm radio flux) of the fitted 
sinusoidals. As one might expect, the quadrupole moment J2 ex- 
hibits no noticeable temporal variation: the sinusoidal fits have 
very low (relative) amplitudes, and the phases inferred from the 
GONG and MDI data are different. There appears to be a sub- 
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Fig. 8. Temporal variation in (even) multipole moments J 2 - Jn- 
The solid lines are the results from using GONG data, dashed 
lines from using MDI data. In all cases a sinusoidal fit, with a 
period of 1 1 years, is also shown. 

stantial difference between the mean value of J2 calculated from 
the GONG and the MDI data. This is most probably a result 
of already known differences in the splitting coefficients (Schou 
et al. 2002). Similar differences between GONG and MDI es- 
timates for J2 were found by Pijpers (1998). We note in pass- 
ing that Emilio et al. (2007), using measurements by MDI, find 
strong variation in the figure of the Sun, which is contrary to 
their earlier results (Kuhn et al. 1998); however, most of the dis- 
tortion from sphericity is the direct response to the centrifugal 
force on the rotating surface layers, and not from the asphericity 
of the gravitational field. The temporal variation in the distor- 
tion observed at the solar surface is certainly associated with the 
magnetic field, which itself is correlated with the zonal bands of 
alternate fast and slow rotation, and is probably confined to lay- 
ers near the surface. Of course, the uncertainty in the estimate of 
the gravitational distortion obtained from the shape of the sur- 
face of the Sun is much greater than the uncertainty with which 
J2 can be obtained from seismic measurements. 

Kuhn et al. (1998) found significant temporal variation in the 
P4(cos 9) component of distortion at the solar surface, although 
their errors are much larger in this component. Accordingly, we 
seek a corresponding variation in J 4. It appears in Figure 8 that 
there is indeed some systematic variation, for inferences from 
both GONG and MDI data behave similarly. However, the am- 
plitude of the sinusoid fitted to the GONG data is only about 
60% of that of the corresponding MDI sinusoid, and, moreover, 
the phases differ by about 1.3 years. The correlation coefficients 
between the inferred values of J '4 and the 10.7 cm radio flux 
are only 0.29 and 0.5 1 for GONG and MDI data respectively, al- 
though these quite low values appear to be a consequence mainly 
of the phase difference. However, the results do suggest that 
J 4 may undergo temporal variation with an amplitude of about 
10" 10 , which is 2% of the mean, and that this variation is cor- 
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Table 1. Temporal variation in multipole moments 



k 




</*> 






Phase (years) 


\SJ, 


■i 








GONG 
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GONG 
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GONG 
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2 


(2.18 4 


: 0.005) x 10~ 7 
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: 0.009) x 10" 7 
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1.0 x 10- 10 


9.5 x 10 
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4.9 x 10-" 
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6 
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+ 0.2) x lO" 10 


-(0.8 


± 0.3) x 10" 10 


7.3 
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4.8 x 10~ 12 
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±0.9)x 10-" 


4.5 


3.6 


1.6 x 10-" 
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1.6 


1.4 


7.3 x 10~ 12 


7.8 x 10" 


-12 


12 


-(1.9 


+ 2.4) x 10~ 12 


-(1.8 


± 2.0) x 10~ 12 


6.8 


6.9 


3.2 x 10~ 12 


2.6 x lO- 


-12 



related with solar activity. The variation is much less than what 
Kuhn et al. (1998) found in the P 4 (cos 8) component of the dis- 
tortion from sphericity of the shape of the Sun's surface; the dif- 
ference could be in part a product of large errors in the distortion 
observed, but it is likely to be due mainly to the direct effect of 
the centrifugal force, and also the magnetic field and accompa- 
nying velocity and temperature variations associated with solar 
activity in the low-density surface layers, which influence the 
figure but hardly modify the gravitational field, even though the 
kernels for the high-order multipole moments are more concen- 
trated near the surface. The mean value of J \ in Table 1 is close to 
that found by Roxburgh (2001) and Mecheri et al. (2004), who 
used rough approximations to the seismically inferred rotation 
rate. 

The variations in 7e that we have inferred from the GONG 
and MDI data do not agree: whereas the MDI results suggest a 
clear variation which is correlated with solar activity, the GONG 
results reveal none. The mean values of J(, differ by a factor 3, 
although they are at least of the same order as that found by 
Roxburgh (2001). 

The variation in Jg is anticorrelated with solar activity, and 
now the variations inferred from both GONG and MDI data look 
similar. However, the mean values do not even have the same 
sign. With Jio and J\2 there is good agreement between GONG 
and MDI. The magnitudes of J% and higher-order moments are 
comparable with their temporal variations, although we hasten to 
add that they are comparable also with the estimated 1-cr errors 
in the individual estimates. Therefore they may not be signifi- 
cant. But we have presented them because they do tend to show 
definite temporal trends which are similar in both the GONG and 
MDI inferences. Higher-order multipole moments also exhibit 
variation, but we do not present them because their reliability is 
certainly questionable. 

We point out that the contribution from rotation to each mul- 
tipole moment can be expressed as a bilinear function of the odd 
degeneracy splitting coefficients of the seismic modes, a func- 
tion which depends also on the density p(r). Therefore, because 
the density of the same reference solar model was used to com- 
pute all the moments, whether the seismic data used were from 
GONG or MDI, all the differences amongst the inferences must 
have their origin solely in the differences amongst the measured 
splitting coefficients. 

5. Summary and Conclusions 

We have studied the temporal variation in the rotational kinetic 
energy and angular momentum of the Sun in different regions of 
the convection zone. We find that in the outer layers (r ^ 0.9/? Q ) 
the variation is well correlated with solar activity. At high lati- 
tudes this correlation continues throughout the convection zone, 



Table 2. Rotational kinetic energy, thermal energy and the esti- 
mated rms magnetic field using Eq. (5) in different layers 



Radial extent 


Mean rot. KE 


rms B 


Thermal energy 




(erg) 


(G) 


(erg) 


0.99 < r/R < 1.00 


5.3 x 10 37 


250 


6.2 x 10 40 


0.95 < r/R Q < 0.99 


6.5 x 10 39 


1400 


2.4 x 10 43 


0.85 < r/R Q < 0.95 


9.3 x 10 40 


3700 


1.4 x 10 45 


0.75 < r/R B < 0.85 


2.0 x 10 41 


6100 


7.6 x 10 45 


0.65 < r/R Q < 0.75 


2.6 x 10 41 


7900 


2.2 x 10 46 



while at low latitudes (and also when integrated over all lati- 
tudes), there is a transition to anticorrelation around r = Q.9R & . 
Through the bulk of the convection zone the high latitudes seem 
to be rotating faster during the solar activity maximum, and the 
equatorial regions rotate more slowly. It is difficult to draw any 
definite conclusions about the temporal variation below the con- 
vection zone, however, because the errors are large. 

The amplitude of the temporal variation in the rotational ki- 
netic energy integrated over the entire convection zone is of or- 
der of 3 x 10 38 erg. This implies a rate of variation whose am- 
plitude is about 5 x 10 30 erg s _1 . It is not clear if this varia- 
tion is compensated for by opposite variation in deeper layers or 
whether it represents exchange of energy between rotation and 
magnetic field or some other form. It is interesting to note that 
this value is roughly comparable to (actually about 3 times) the 
amplitude of the observed variation in the solar irradiance, which 
is about 0.04%, although the correspondence might well be for- 
tuitous. It is, however, tempting to speculate that the temporal 
variation in rotational kinetic energy, and the magnetic energy, 
has some role to play in the irradiance variation. A part of the 
variation in the rotational kinetic energy would be expected to 
be compensated by the magnetic energy stored in the solar con- 
vection zone. If that were the case at all radii, we can estimate 
the spatial variation of the rms magnetic field using the relation 

d2 

— 4nr 2 5r * ST. (5) 
8tt 

The total variation 5T in T over the solar cycle is about 0.002 T, 
where T(r\ , ri) is the mean rotational kinetic energy in the region 
between radii r\ and r^. Using this, together with the mean ki- 
netic energy integrated over all latitudes within several spherical 
shells, we list in Table 2 these admittedly very crude estimates 
of the magnetic field. 

The gravitational quadrupole moment J2 does not exhibit 
any significant temporal variation, whereas some other, higher- 
order, multipole moments do. The variation in J4 is correlated 
with solar activity, but the amplitudes of variation seen in GONG 
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and MDI data are somewhat different. Clear variation in 7g - Jn 
is inferred from both GONG and MDI data, the variation in J\q 
being correlated with solar activity and that in J$ and J\2 anti- 
correlated. The absence of temporal variation in J2 is consistent 
with observed absence of variation in the P2(cos 6) component 
of distortion at the solar surface (Kuhn et al. 1998). This is to be 
expected, since J2 is determined essentially by the spherically 
averaged component of the angular velocity, which suffers lit- 
tle relative variation. The higher-order multipole moments are 
determined also by the latitudinally varying components of the 
angular velocity, which are known to suffer temporal variation 
correlated with solar activity. Note that the predominant con- 
tribution to J2 arises from regions near the base of the convec- 
tion zone (Gough, 1981; Pijpers 1998), whereas the higher-order 
multipole moments arise largely from contributions from outer 
layers of the Sun, where the temporal variation in the rotation 
rate is not insubstantial. Of course, the magnitude of J2 is much 
larger than that of the higher multipole moments, and it is likely 
that J2 also has temporal variations with comparable absolute 
magnitude. 
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